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ABSTRACT 

Many open clusters have a deficit of observed white dwarfs (WDs) compared with predictions of the 
number of stars to have evolved into WDs. We evaluate the number of WDs produced in open clusters 
and the number of those WDS detectable using photometric selection techniques. This calculation 
includes the effects of varying the initial- mass function (IMF), the maximum progenitor masses of 
WDs, and the binary fraction. Differences between the calculated number of observable WDs and 
the actual number of WDs observed in a specific cluster then indicate the true deficit of WDs that 
must be explained through effects such as dynamical evolution of the cluster or close binary evolution. 
Observations of WDs in three open clusters, the Hyades, Pleiades, and Praesepe, are compared to the 
calculated observable populations in those clusters. The results suggest that a large portion of the 
white dwarf deficit may be explained by the presence of WDs in unresolved binary systems. However, 
the calculated WD populations still over-predict the number of observable WDs in each cluster. 
While these calculations cannot determine the cause of this residual white dwarf deficit, potential 
explanations include a steep high-mass IMF, dynamical evolution of the cluster, or an increased 
likelihood of equal-mass components in a binary system. Observations of complete WD samples in 
open clusters covering a range of ages and mass can help to distinguish between these possibilities. 
Subject headings: white dwarfs — open clusters and associations: general — open clusters and asso- 
ciations: individual (Hyades, Praesepe, Pleiades) — binaries: general 



1. INTRODUCTION 

White dwarfs (WDs) are the final endpoint of stellar 
evolution for the vast majority of stars. WDs are gaining 
importance in tracing the history of stellar populations in 
the g alaxy, including the galact ic disk (e.g.jWinget et alJ 
Il987|) . open star cluste rs (e.g. IRicher et alJll998j) . and 
globular star clusters ((Hansen et alJ 2002). WDs are 
also useful in studies of supernova physics, as the upper 
mass of WD progenitors represents the critical progenitor 
mass (Merit) 1 for core-collapse supernovae. The value of 
M cr it is relatively uncertain, with the best estimates of 
5.5M Q < M crit < 10M Q . 

The first significant population of white dwarfs (WDs) 
in ope n clusters was identified by lEggen fc Greensteinl 
(1965) i n the Hyades, Praesepe, and the Pleiades. 
iTinslevI l|1974|) compared the number of Hyades WDs to 
the estimated number of Hyades stars having completed 
their lifetimes. Her work found that the WD numbers 
are compatible with M cr it ~ 4Mq, though this determi- 
nation depends sensitively on the assumed initial-mass 
function (IMF), the turnoff mass of the Hyades, and the 
Hyades distance modulus. This value of M cr it is lower 
than current lower limits [Merit ^ 6M©), set in large 
part by the existence of LB 1497, the lone WD in the 
P leiades (turnoff mass ~ 5 .4M^). 

iWeidemann et al.1 l|1992|) discuss the Hyades WD pop- 
ulation and determine a deficit of 21 WDs (28 WDs pre- 
dicted versus 7 observed WDs) for M cr it — 8M Q and 
ascribe this white dwarf deficit to dynamical evapora- 
tion of WDs from the Hyades. Similar deficits of WDs 
have since been observed in other open clusters, including 
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M67 IRicher et al.lll998l) . Praesep e (Claver et al 2001 , 
hereafter ICLBKD . and NGC 2099 ijKalirai et alJl2001c|) . 
The existence o f the white dwarf d eficit is still debatable, 
as discussed bv lvon Hippell iri998|) . 

Three explanations for the white dwarf deficit in 
open clusters readily come to mind. First, WDs 
may evaporate from open clusters due to dynamical 
evolution. Mass segregation and galactic tidal fields 
alone are not sufficient to remove WDs from open 
clusters, as a 0.6M© WD is still more massive than 
typical cluster stars and thus unlikely to suffer pref- 
erential evaporation. This is borne out by modern 
jV -body simulations of open cl usters, such as those 
of iPortegies Zwart et al.l (120011): iBaumgardt fc Makinol 
(I2003|) and lHurlev Sharal (|2003Il which find that WDs 
remain bound in open clusters. However, if a WD re- 
ceives a sufficient velocity kick from asymmetric mass 
loss during its post-main sequence evolution, the WD 
may become unbound from the open cluster. This sce- 
nario was first suggested bv IWeidemann et al.l l|1992|) to 
explain the WD defi cit in the Hyades . Rece nt, simple N- 
body simulations by Fellha uer et al.l 1)200.11) confirm that 
this mechanism can preferentially remove WDs from an 
open cluster, though scant observational evidence exists. 

A second explanation for the white dwarf deficit, 
though not exclusive of the dynamical evolutionary argu- 
ment, is that the WDs may be hidden in binary systems 
in which the intrinsically faint WDs are not detected 
due to the overwhelmin g light of brighter companion 
(e.g. iKalirai et af1l2001c() . Searches for WDs in binaries 
have been undertaken in the Hyades. Two of the known 
Hyades WDs, EGGR 38 and V471 Tau, are in binary 
systems. More recently, additional Hyades WDs have 
been discovered hidden in unresolved binaries, includ- 
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ingHD 27483 (|Bohm-Vitenselll993lL VA351 <|Franz et all 

1998), and four pot ential WDs in Am binary star systems 
Ipebernardi et aQ[200C| . Clearly some WDs lie hidden 
in unresolved binary systems, but the exact numbers are 
not known. 

The third explanation for the small number of white 
dwarfs seen in the Hyades and other well studied young 
clusters could be that our expectations of WD numbers 
are incorrect. The deficit discussed for the Hyades would 
go away if M cr n is low (~ 4M ), if the IMF for masses 
above the cluster turnoff mass is steeper than the present- 
day mass function around the turnoff mass, or a combi- 
nation of these two effects. The low value of M cr n seems 
unlikely, given that WDs have been observed in clus- 
ters with turnoff ma sses higher than 4M , inc luding the 
Pleia des. NGC 2516 (IKoester fc Reimersll996j) and NGC 
2168 l|Reimers fc Koesterlll988ft . 

We have developed a Monte Carlo method of calcu- 
lating the number of WDs detectable in observations of 
specific open clusters, a calculation designed to aid in 
the interpretation of data from the Lick- Arizona White 
Dwarf Survey (LAWDS, iWilliams fc BoTtel 12003]) and 
could be used in conjunction with other ongoing clus- 
ter WD surveys, such as the CFHT Open Cluster Sur- 
vey (iKalirai et alJl2001al) and the WI YN Open Cluster 
Study l)von Hipnel fc Saraie dini 1998). This calculation 
utilizes the observed characteristics of specific open clus- 
ters, including the cluster age, distance and reddening, 
and determines how many WDs would be detected given 
photometric WD selection criteria. The structure of this 
calculation, its limitations, and tests of the calculation 
are presented in §2. In §3 we present the results of the 
calculations for the Hyades, Praesepe, and the Pleiades 
and compare these results to the observed WD popula- 
tions. In §4 we discuss the results and discuss the use- 
fulness of these calculations in regard to current searches 
for WDs in open clusters. 

2. CALCULATING THE WD POPULATIONS IN OPEN 
CLUSTERS 

The number of WDs detectable by photometric obser- 
vations of open clusters can be calculated by appropriate 
modeling of the cluster, including the effects of unre- 
solved binary stars, choice of IMF, WD detection criteria, 
M cr it, metallicity and cluster age. Assuming that such 
an open cluster model is sufficiently realistic, any signif- 
icant difference between the calculated numbers of ob- 
servable WDs and the actual numbers of observed WDs 
represents the "true" white dwarf deficit that must be 
explained by other means, such as evaporation or other 
dynamical evolution. 

Monte Carlo techniques are used to calculate the ob- 
servable WD population. Stars are randomly drawn from 
an input IMF. For this work, four input IMFs £(M) oc 
M~( 1+r ) were c onsidered: a Salpeter IMF with T — 1.35 
l)Salpeterl lT955j). a steeper power law that accounts for 
unres olved binary systems with T = 2 l)Navlor et al.l 
2002), and a broken power law with a flat low-mass 
slope: T = 0.2 (M < 1M Q ) ; T = 1.8 (M > 1M ) 
ijNavlor et alJl2002|) . hereafter called the "Naylor IMF," 
and the Kroupa IMF with T = 0.3 (M < 0.5M Q ) ; T = 
1.3 (A/ > O.5M ) IKrounal ll200T1) . 

Stars are assigned a binary companion with a prob- 
ability based on an input binary fraction; the binary 



mass fraction is assumed to be random, i.e. the pri- 
mary and secondary stars are drawn from the same IMF 
(see fl2.ll below). The stars are then evolved to the 
input cluster age using the Z = 0.008, Z = 0.019 or 
Z = 0.04 solar-scaled metallic ity stellar evolutionary 
models of Sal asnich et alJ ((2000). These models provide 
evolutionary data on stars with masses up to 2OM , suf- 
ficiently massive to explore all likely values of M cr it- 

Two outcomes are assumed for stars having completed 
their evolution. Stars with > M cr it undergo super- 
novae. Stars with M* < M cr u form a WD. The WD 
mass is determined from a linear approximati on of the 
initial- final mass relation in iWeidemannl l)2000D . 

M W d = |r M * + O.414M . (1) 
35 

Each WD is assigned a cooling age equivalent to the pro- 
genitor lifetime subtracted from the input cluster age. 
The radii, effective temperature, and surface gravities of 
each WD are interpolated fr om th e carbon/oxygen mod- 
els of Wood (1995, hereafter IW95l) . 

Photometric indices for each system are determined 
by summation of the flux from each component. Super- 
nova remnants are assumed to emit no flux. WD fluxes 
are determined from the updat ed model atmospheres of 
iFinlev. Koester. fc Basril l)1997l) graciously provided by 
D. Koester. Stel l ar flux es are ta ken from the iso c hrones 
of iCirardi et all |2Q02) for the ISalasnich et"afl (2000) 
evolutionary models. These indices are corrected for the 
input distance modulus a nd reddening, assuming th e in- 
terstellar reddening law of lRieke fc Lebofskvl l)1985|) and 
Rv = 3.1. 

Once the photometric indices for each system have 
been calculated, these indices are compared to the in- 
put WD photometric-selection criteria. These criteria 
include color indices and single-band magnitude limits, 
with the intention that the photometric criteria used to 
select WD candidates in actual observations of a given 
cluster are used in the simulation of that cluster. If a 
system meets these input photometric selection criteria, 
the number of detected WDs is incremented by one. If it 
fails to meet the selection criteria, the WD is considered 
undetected. 

The Monte Carlo calculation continues to select stel- 
lar systems until the number of stars brighter than a 
limiting magnitude Vu m reaches the input value. For ex- 
ample, if the cumulative luminosity function of an actual 
open cluster has AT* stars brighter than Vn m , then each 
simulated cluster will also contain AT, stars brighter than 
Vum- 10 4 realization of each cluster are used to ensure 
high precision in the output. 

2.1. Limitations of the calculations 

These calculations represent a simplified version of ac- 
tual open clusters. The simulation assumes that all bi- 
nary stars are unresolved, an assumption whose validity 
depends strongly on the cluster distance and the distri- 
bution of physical separations among binary stars. Bi- 
nary interactions during stellar evolution, such as mass 
transfer systems, common envelope stages, cataclysmic 
variables, and the potential disruption of binaries due 
to mass loss are ignored. Open cluster dynamical evo- 
lution, such as evaporation of stars and evolution of the 
binary fraction, is also ignored. The binary mass ra- 
tios are assumed to be random, which is likely not be 
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the case for close binaries, but is reasonable for wide r 
binaries {Mason et aLllToMIAbt. Gomez, fc Levvl [l990). 
Finally, stellar systems with more than two components 
are not considered. 

Ignoring binary-star evolution and the different mass 
ratios of close binaries than a random IMF results 
in far fewer double-degenerate systems being created 
than as seen in realistic iV-body simulations, such as 
the significant double -degenerate population seen in 
IHurlev fc Sharal lEffll . The presence of any observed 
double-degenerate sequence in an open cluster therefore 
can provide a constraint on the relative importance of 
close-binary evolution in open clusters. 

The lack of consideration of dynamical evolution is not 
a deficiency in the calculations, but rather is one purpose 
of these calculations - to determine the actual number of 
missing WDs that dynamical evolution must be invoked 
to explain. The strong, observed evaporation of low-mass 
stars has no great impact on this calculation. Low-mass 
stars in binaries with a (comparatively) high-mass WD 
should not be affected greatly by evaporation processes, 
and low-mass stars not in a binary system with a WD are 
not used in the calculation beyond the initial selection of 
cluster stars. 

In short, these calculations are not a simulation of WDs 
in the larger sense of open cluster evolution. Rather, 
this calculation provides a more sophisticated means of 
comparing the observed WD populations of specific open 
clusters to the expected observable WD population. Any 
differences between the actual observed WD population 
and the calculated observable population provide a mea- 
sure of the strength of the dynamical evolution of the 
WD population. 

2.2. Trials of the calculations 

We ran numerous tests to verify that these calculations 
return reasonable results for certain limiting cases. Fig.^ 
shows the output IMF of a single simulated cluster with 
5000 stars more massive than 1M© for each of the input 
IMFs. As one would hope, the output IMFs agree with 
the input IMFs. 

Simulated clusters with a T — 2 power-law or the Nay- 
lor IMF are found to contain fewer WDs than those with 
a Salpeter or Kroupa IMF, merely due to the fact that 
fewer high-mass stars formed due to the steeper slope. 
The Naylor and Kroupa IMFs result in a higher fraction 
of WDs being hidden in binary systems than the single 
power-law IMFs due to the relatively flat low-mass IMF 
slope, as a companion to a WD is more likely to be higher 
mass and therefore more luminous than the WD. 

Fig.|21shows output color-magnitude diagrams (CMDs) 
for calculations of clusters at four different ages contain- 
ing no binary stars. Occasional minor numerical issues 
arise during the interpolation of photometric indices for 
evolved stars due to the large changes in magnitudes with 
minute changes in stellar mass. These issues arise pri- 
marily at inflection points in the isochrones, such at the 
helium flash and at the tip of the AGB. These errors 
do not affect the simulated detection of WDs, as both 
the correct and erroneous magnitudes for these stars are 
bright enough to hide any WD companions, and the num- 
ber of stars affected is very small. 

The effect of the addition of binary stars is shown in 
Fig. |2| The sequence of binaries consisting of two main- 
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Fig. 1. — Simulated initial-mass functions. Histograms indi- 
cate the output for a Salpeter IMF (solid), a T = 2 power-law 
IMF (dotted), Naylor IMF (long-dashed) and Kroupa IMF (dash- 
dotted). Lines indicate the expected number of stars for each IMF. 
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Fig. 2. — B — V, My CMD for simulated clusters of various ages. 
Simulated ages are for logarithmic cluster ages of 7.5 (pluses), 8.0 
(asterisks), 8.5 (open circles), and 9.0 (filled squares). Simulated 
clusters have 1000 existing stars more massive than IMq. 



sequence (MS) stars is readily visible in the B — V, My 
CMD. The fan of points leading from the MS to the WD 
cooling sequence is composed of MS- WD binaries. The 
B — V,U — B color-color diagram illustrates the advan- 
tage of three-band photometry over two-band photom- 
etry. The WD-MS binaries, which would blend in with 
background disk MS stars in a B — V, V CMD, are well- 
separated from the locus of main-sequence in color-color 
space. In addition, the color-color space information 
clearly separates the WD cooling sequence from massive 
MS stars. 

A WD- WD binary sequence is visible only when the 
number of stars in a cluster is very high, as seen in the 
lower-right panel of CMD in Fig. |3| Any significant ob- 
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Fig. 3. — B — V , My CMD (left) and B — V, U — B color-color diagram (right) for simulated clusters with binary stars. For each diagram, 
the number of stars more massive than IMq is 1000 (top panels), 10 4 (bottom left) and 10 5 (bottom right). The top left panel in each 
diagram is for a simulated cluster with no binary stars; the binary fraction for each of the other panels is 0.5. The WD-WD binary sequence 
is visible in the lower right panel of the CMD, 0.75 mags brighter than the lone WD cooling sequence. Each simulation is of a cluster with 
a logarithmic age of 8.50 and the broken power-law IMF. 



served WD-WD binary population in open clusters would 
therefore confirm that the binary component mass ratios 
are not random and/or that dynamical processes such as 
binary interactions are not negligible. iHurlev fc~ Shara 
( 2003) find evidence of a WD-WD sequence i n the pho- 
tometry of NGC 6819 of lKalirai et"aiT l|2001bj) . 

Multiple (10 4 ) realizations of each cluster are used to 
calculate precise means and standard deviations for var- 
ious categories of WDs, including the total number of 
WDs, the total observed number of WDs, the number of 
WDs in binary systems, etc. 

3. SIMULATIONS OF SPECIFIC CLUSTERS 

We have calculated the predicted observed WD popu- 
lations for three open clusters: the Hyades, the Praesepe, 
and the Pleiades. The output populations are compared 
to the number of observed WDs in each system. The 
calculated numbers of WDs in each cluster fulfilling var- 
ious criteria, including whether the WD met the photo- 
metric selection criteria (these are given individually for 
each cluster below), whether a WD was in a binary sys- 
tem, and WD progenitor masses were calculated, and the 
mean and standard deviation of each over the ensemble 
of trials were calculated. 

The standard deviations are found to be nearly equal to 
those expected from Poisson statistics. Therefore, we use 
Poisson statistics to compare the calculated populations 
to observations of the actual cluster WDs. The probabil- 
ity P of obtaining the actual number of observed WDs in 
a category or fewer is calculated from Poisson statistics 
given the mean value obtained from the Monte Carlo cal- 
culation. If P < 0.05, then the derived WD population 
is considered to be inconsistent with the observed popu- 
lation. For probabilities 0.15 > P > 0.05, the calculated 
is termed to be mildly inconsistent with the data. 

The Hyades — The Hyades is the open cluster with the 
best-studied WD system. Nine WDs are traditionally 



considered cluster members, though some questions as 
to the membership of s ome of these have been raised 
ijWeidemann et alJH992|) . Of these, the seven lone WDs 
and the WD-MS binary EGGR 38 are detectable as blue 
excess objects in t/Byphotometry. 

For the Monte Carlo realizations the Hyades, the fol- 
lowing physical parameters are used. The intrinsic dis- 
tance modulus is assumed to be 3.33 llPerrvman et alJ 
1998), the color excess E{B-V) = 0.01 ijMermTmodl 
1995), and the logarit hmic age is assumed to be 8.80 
ijPerrvman et alJll998[) . The binary star fraction is cho- 
sen to be 0.4, approximately the lower limit of 0.41 ±0.05 
determined from the combinatio n of speckle, spectro - 
scopic, and direct imaging surveys l|Patience et al.ll 998). 

The Hyades WDs have been detected and confirmed 
through a variety of photometric, spectroscopic, and as- 
trometric surveys, and a direct comparison of these se- 
lection criteria with the calculated population is not per- 
formed. In order to determine if a cluster WD is to be 
considered observable, we apply the photometric criteria 
of U-B < 0.0 and B — V < 0.6 to select WDs. These 
values encom pass the range of indices o bserved in the 
WD sample of lEggen fc Greensteinl l) 19651) , including the 
seven single Hyades WDs and the binary EGGR 38. The 
MS- WD binaries not selected via these criteria represent 
the contribution of binaries to the WD deficit. 

The metallicity of th e Hyades is super-sol ar, with 
[Fe/H] = 0.14 ± 0.05 ijPerrvman et all 11998(1 . corre- 
sponding to Z = 0.034 for a solar Z = 0.0245 
l)Christensen-Dalsgaard et al]ll996|K We therefore cal- 
culate WD populations for each of two mctallicitics, 
Z = 0.019 and Z = 0.040, for comparison. The cal- 
culated population also explores three values of M cr a, 
10M Q , 8Mq, and 6M Q , bracketing the range of likely 
values. Finally, the calculations were performed using 
each of the four IMFs discussed above. 

Each realization of the cluster contains 138 stellar sys- 
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terns in the magnitude range 3.0 < B < 10.0, c onsist ent 
with the Hyades luminosity function of lOortl l|1979l) . a 
complete luminosity function of the Hyades upper main 
sequence. The resulting WD populations of the Hyades 
are given in Tabled 

Our results indicate that the eight photometrically- 
selected Hyades WDs are inconsistent (P < 0.03) with 
what would be expected from a Salpeter or Kroupa IMF, 
independent of the value of M crit . However, the calcu- 
lated number of WD detections for the r = 2 power- 
law IMF and the Naylor IMF are fully consistent with 
the observed number of WDs for Z = 0.019, and con- 
sistent (P « 0.2) for Z = 0.040. The lack of any 
observed Hyades WDs with pro genitor masses > 4M 
ijWeidemann et aT] 1 19921 ICLBK1 is mildly inconsistent 
for the r = 2 power law and Naylor IMFs at M cr i t = 
6M Q (P = 0.13 and P = 0.07, respectively) and for the 
r = 2 power-law IMF for M crlt = 8M© (P = 0.07), but 
is inconsistent with the survey results for the remaining 
explored parameter space. This result will be explored 
in the discussion (§4) below. 

Of the eight photometrically-detectable Hyades WDs, 
one is in a binary system. This is inconsistent with the 
number of binaries detectable from a Salpeter IMF (P < 
0.05) and mildly inconsistent with the Kroupa IMF and 
r = 2 power-law IMF for the binary fraction of 0.5, but 
is otherwise consistent with the calculations for the Z = 
0.019 systems. For Z = 0.04, the single observed binary 
is inconsistent with one observed binary for each IMF 
except the Naylor IMF, which is fully consistent with 
t he observations 

Boh m-Vitensd l)1995l) reports on an IUE survey for WD 
companions to F stars in the Hyades, with only one WD- 
F star system (actually a triple system with two F stars 
and a WD) being detected. In our calculations, we have 
identified all stars meeting the B — V selection criteria in 
that study. The average number of WD companions to 
these F stars are presented in the final column of Table 
On average, < 1 F star- WD binary is present in a 
simula ted cluster, fully consistent with the lBohm-Vitensel 
lfl995h results. 

Praesepe — ICLBKI list ten candidate Praesepe WDs. 
Based on their analysis of six spectra, five are cluster 
WDs, and the other is a confirmed non-member. We ob- 
tained spectroscopy of two of the remaining candidates, 
WD 0834+209 and LB 1839, w ith the upgraded blue 
camera of LRIS llOke et al.lll995ft on the Keck 10-m tele- 
scope in December 2002. Both of these objects are found 
to be QSOs, with LB 1839 at z = 0.86 and WD 0834+209 
at z = 1.09. Therefore, at most seven of the Praesepe 
candidates are WDs, and this is the number of Praesepe 
WDs used in the comparison to the calculated popula- 
tions. 

The realizations of the Praesepe calc ulations were nor- 
malize d by the luminosity function of I.Tones fc Staufferl 
(1991), who find 179 stars (com pleteness-correcte d ) with 
2.5 < My < 7.5. While the I.Tones fc Staufferl |l99lT) 
data extend to fainter magnitudes, CLBK find that the 
Padova isochrones do not match the Praesepe MS below 
this magnitude, and this mismatch would systematically 
bias our calculations. 

The age and metallicity of the Praesep e are ne arly in- 
distinguishable from that of the Hyades ( CLBK) , so we 
adopt the Hyades logarithmic cluster age of 8.80 for Prae- 



sepe. As the metallicity was found to have only a small 
effect on the calculated observable number of WDs in the 
Hyades, we simulate Praesepe using only the Z— 0.019 
isochr ones. The binary star fract ion of the Praesepe is 
~ 0.4 l|Mermilliod fc Mavorlll999l) . so this value is used 
in the simulation. 

For WD detection criteria, we requir e U — B < 0.0 and 
B — V< 0.6, limits encompassing the ICLBKI WD sam- 
ple. As in the Hyades, the form of the IMF and M ar it 
were varied. The calculated WD populations are given 
in Table H 

As with the Hyades, the observed WD numbers are 
inconsistent with the calculated output for a Salpeter 
and Kroupa IMFs. The numbers of observable WDs cal- 
culated using the Naylor law IMF are consistent with 
the observed number for M cr u — 6Mq (P = 0.24) but 
mildly inconsistent for higher values of M cr i t (P = 0.15 
and P = 0.11 for M crit = 8M Q and M crit = 10M o , re- 
spectively) . The r = 2 power-law IMF simulation results 
a re fully consistent with the observed numbers. 

CLBK] detected one WD with a progenitor mass > 
4Mq. This number is inconsistent with a Salpeter and 
Kroupa IMF for all values of M crit (P < 0.04). The 
simulated number of WDs with massive progenitors is 
inconsistent with the Naylor IMF for M crit = 10 M Q 
(P = 0.05), but consistent with the one observed WD 
for M crit = 8M Q and M crit = 6M© (P = 0.17 and 
P = 0.20, respectively). The r = 2 power law WD num- 
bers are reasonably consistent with the observations for 
all values of M cr i t . 

None of the five confirmed Praesepe WDs are known 
to be in binary systems. If neither of the two remaining 
candidates is in a binary system, then the lack of binaries 
is inconsistent with the simulated results for a Salpeter 
or Kroupa IMF. The lack of binaries is mildly inconsis- 
tent (P = 0.08) with the T = 2 IMF, and consistent 
(P w 0.19) with the Naylor IMF. It is emphasized that 
this inconsistency is between the lack of observed WD 
binaries in Praesepe and the number of simulated WDs 
in binaries which would be detected given the stated pho- 
tometric selection criteria; this discrepancy does not in- 
clude the WDs in binary systems with photometric in- 
dices outside the selection criteria. 

PleAades — The Pleiades has one known member WD, 
LB 1497. For this cluster, we have normalized the 
IMF such that each simulated cluster contains 107 stars 
in the magnitude range 2.0 < V < 10.0, con sistent 
with the luminosity function of the Pleiades from [ Jones! 
(QUI. We assume the iPinsonneault et all l)1998j) dis- 
tance modulus of 5.60 and reddening of E(B — V) — 
0.04. The metallicity of the Ple iades is near so- 
lar, w ith [Fe/H] = -0.034 ± 0.024 l|Boesgaard fc Friel 
11990ft. so we use the Z=0.019 isochrones. We adopt 
the iMermilliod et alJ Jl992) multiple-star frequency of 
0.36 as the binary frequency. Age estimates of the 
Pleiades have s hown scatter, with recent estimates be- 
tween 100 Mvr IPinsonneault et all 1 19981) and 125 Myr 
l|Stauffer. Schultz. fc Kirkpatrickl I1998D . We therefore 
produce calculations for logarithmic cluster ages of 8.0 
and 8.1. 

The results are given in Table OH For M cr u = 8M 
and M cr it = 10-M©, the Salpeter and Kroupa IMFs are 
inconsistent with the one observed WD, regardless of the 
adopted cluster age. For M cr it — 6M Q , the Salpeter and 
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Kroupa IMFs are mildly inconsistent with the observed 
WD if the Pleiades have a logarithmic age of 8.1 (P = 
0.12), but are consistent for a logarithmic age of 8.0. The 
T = 2 and Naylor IMFs produce WD populations mildly 
inconsistent with the Pleiades for M cr u = 10M Q and for 
the M cr it = 8M0, log(age) = 8.1 simulations. Otherwise, 
these two IMFs result in WD populations consistent with 
the single observed WD. 

4. DISCUSSION 

A few observations can be made from the compar- 
isons between the calculated WD populations and the 
observed populations presented above. The calculations 
consistently over-predict the number of observed, color- 
selected WDs and the number of color-selected WDs in 
binary systems. However, the discrepancy between the 
expected number of WDs and the observed number is 
lower than in previous observational studies due to the 
inclusion of binaries in the simulation. Therefore, pref- 
erential evaporation of WDs from these clusters due to 
dy namical evolution need n ot be as severe as suggested 
bv lWeidemann et alJ (^)92) but cannot be ruled out. It 
should also be noted that the binary fractions used in this 
work represent the high end of the range of commonly- 
quoted cluster binary populations. A lower actual binary 
fraction would result in a larger WD deficit. 

The number of observed WDs can be brought into 
agreement with the calculated numbers via at least four 
different methods. First, the slope of the IMF could be 
steeper than the slopes used here. Second, the binary 
fractions may be higher, resulting in more WDs being 
hidden in unresolved binaries. Third, the binary mass 
ratio could be closer to unity than the random parings 
discussed here, which would result in more binary WDs 
being hidden in systems. Fourth, some sort of dynamical 
evolution may be removing WDs from the open clusters. 
From the simulation alone, it is not possible to differen- 
tiate between these scenarios. 

As mentioned above, the lack of observed WDs in the 
Hyades and Praesepe with progenitor masses > 4M Q is 
inconsistent with the calculated populations. A steeper 
high-mass IMF would explain the lack of these WDs, 
as would dynamical evolution. However, there is some 
evidence that binary stars in this mass range tend to 
have a fairly flat distributio n of mass ratios, especially 
for orbital periods < O.lyr (jAbt. Gomez, fc Levvlll9961: 
iMason et al.l ir9981. If this is true, WDs with massive 
progenitors may be more likely to be in binary systems 
with massive, unevolved stars than in binary systems 
with fainter, low-mass stars, and therefore more likely to 
go undetected. 

There are several avenues of research which could shed 
some light on these issues. First, systematic, thorough 
searches for WD-mass companions to stars in the Hyades 
and Praesepe are needed to complete the census of WDs 



in each cluster, thereby determining whether or not bi- 
nary systems can account for the majority of "missing" 
WDs in open clusters. Such a survey could include a 
combination of spectroscopic and high-resolution imag- 
ing (e.g. speckle and adaptive o ptics surveys) o f these 
clusters, such as those published in lPatience et al.l (1998) 
and similar surveys. These surveys could also determine 
if the binary mass ratios for any close binaries are non- 
random, as opposed to the random ratios assumed in this 
work. 

Also, WD searches in open clusters of a wide range of 
ages, such as the LAWDS survey we are currently un- 
dertaking, can provide evidence as to whether the deficit 
of massive WDs is increasing with time, as would be 
expected for dynamical evolution. In addition, detailed 
comparisons of clusters of similar ages but differing com- 
pactness and richness should also show differences in the 
white dwarf deficit if dynamical evolution of the WD 
population is significant. 

The calculations presented here will be useful in inter- 
preting data from the WD observational programs. Most 
importantly, the calculations permit us to estimate how 
strong the white dwarf deficit is in a cluster, given basic 
assumptions about the IMF and binary fraction. This 
provides a better estimate on the significance of any ap- 
parent deficit and the cluster-to-cluster variations in the 
deficit than could be gleaned from simpler WD popula- 
tion estimators, such as simple integration of the IMF. 
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Table 1. Simulation of the Hyades. 





z 


Binary 


IMF a 




Nwd 


Nwd 


Nwd 


N W D,F b 


(Me) 




Fraction 








in binaries 


(Mi > 4M Q ) 


g 


0.019 


0.4 


s 


total 


29.7 


16.9 


14.2 


0.4 












17.5 


4.7 


6.7 










p 


total 


14.1 


8.0 


5.6 


0.07 












9.2 


3.2 


2.7 










N 


total 


1 8 Q 
icy 


10 7 


O.U 


U.o 












9.8 


1.7 


3.6 










K 


total 


31.9 


18.0 


15.4 


0.8 










UUoL.l vuu 


16.9 


3.3 


7.1 








0.5 


s 


total 


29.9 


19.7 


14.2 


0.5 












15.5 


5.5 


5.4 










p 


total 


14.2 


9.5 


5.6 


0.08 










U USLl V L.U 


8.4 


3.7 


2.2 










N 


total 


19.2 


12.6 


8.1 


0.9 










u uau v lli 


8.4 


2.0 


2.9 










K 


total 


32.0 


21.1 


15.5 


0.9 










L> UoL.l V LAJ. 


14.5 


3.8 


5.7 






0.040 


0.4 


g 


total 


32.2 


18.2 


15.1 


0.5 










L) UbL.I V LIT 


20.7 


7.0 


7.6 










p 


total 


15.7 


9.0 


6.2 


0.077 










o riQ( 1 r^/oH 

U L/SL^l V LU 


11.3 


4.6 


3.2 










N 


total 


20.6 


11.6 


8.6 


0.8 










UUoL.l V Lll 


11.3 


2.5 


4.0 










K 


total 


34.4 


19.3 


16.4 


0.9 










KJ UbL.I V L.t 1 


19.6 


4.9 


7.9 




U 


n mo 


n A 


o 
o 


LiOLrdl 


25 4 


14 4 




n zl 










(J USL.l V L.t 1 


15.4 


4.5 


4.7 










p 


total 


12.7 


7.2 


4.2 


0.064 










UU5L.1 V L^Ll 


8.5 


3.1 


2.0 










N 


total 


16.8 


9.4 


5.8 


0.7 










observed 


8.9 


1.7 


2.7 










K 


total 


27.1 


15.3 


10.6 


0.7 










observed 


14.7 


3.1 


5.0 




10 


0.019 


0.4 


S 


total 


32.2 


18.2 


16.6 


0.5 










observed 


18.6 


4.8 


7.7 










P 


total 


14.9 


8.5 


6.3 


0.076 










observed 


9.6 


3.2 


3.0 










N 


total 


20.0 


11.3 


9.1 


0.87 










observed 


10.3 


1.7 


4.1 










K 


total 


34.5 


19.4 


18.1 


0.8 










observed 


18.1 


3.3 


8.3 





a S=Salpeter, F = 1.35; P=power law, F = 2; N=Naylor IMF, T = 0.2(M < 1M ) , T = 1.8(M > 1M©); K=Kroupa IMF.r = 0.3(0.080 > 
M > 0.5M Q ), T = 1.3(M > 0.5M Q ) 
b Number of WDs in binaries with F stars (see text) 



Table 2. Simulation of the Praesepe. 



M cr it 


IMF a 


Nwd 


N W d 


N W d 


N W d 


N W d 


Nwd 


(Me) 




total 


observed 


total 


observed 


total 


observed 








in binaries 


in binaries 


(Mi > 4M ) 


(Mi > 4M ) 


10 


S 


35.1 


19.7 


19.8 


4.8 


18.0 


8.4 




P 


13.2 


8.1 


7.6 


2.6 


5.6 


2.7 




N 


23.0 


11.5 


12.9 


1.7 


10.4 


4.7 




K 


38.7 


20.3 


21.8 


3.7 


20.2 


9.3 


8 


S 


32.5 


18.5 


18.4 


4.7 


15.5 


7.3 




P 


12.6 


7.9 


7.2 


2.6 


5.0 


2.4 




N 


21.8 


10.9 


12.3 


1.7 


9.3 


3.2 




K 


35.7 


18.9 


20.1 


3.6 


17.3 


8.0 


6 


S 


27.8 


16.3 


15.8 


4.5 


10.7 


5.1 




P 


11.3 


7.3 


6.5 


2.5 


3.7 


1.9 




N 


19.2 


9.8 


10.8 


1.6 


6.7 


3.0 




K 


30.3 


16.5 


17.1 


3.5 


11.9 


5.6 



a S=Salpctcr, F = 1.35; P=power law, T = 2; N=Naylor IMF, F = 0.2(M < 1M ) , F = 1.8(M > 1M ); K=Kroupa IMF,r = 0.3(0.080 > 
M > O.5M ), T = 1.3(M > O.5M ) 



White dwarfs in binaries 



Table 3. Simulation of the Pleiades. 





log(age) 


IMF a 


%D 


N W D 


N W D 


N W D 








total 


observed 


in binaries 


observed in binaries 


10 


8.0 


s 


9.1 


7.6 


4.8 


3.3 






p 


3.9 


3.5 


2.1 


1.7 






N 


5.2 


3.6 


2.7 


1.1 






K 


9.8 


7.2 


5.1 


2.6 




8.1 


g 


11.3 


9.1 


6.0 


3.9 






p 


4.9 


4.4 


2.6 


2.0 






N 


6.6 


4.4 


3.5 


1.3 






K 


12.1 


8.6 


6.4 


2.9 


8 


8.0 


s 


6.7 


5.7 


3.6 


2.5 






p 


9 Q 


9 7 


1 fi 


1 .o 






N 


3.9 


2.8 


2.0 


0.9 






K 


7.2 


5.3 


3.8 


2.0 




8.1 


S 


8.9 


7.3 


4.7 


3.2 






p 


4.1 


3.6 


2.2 


1.8 






N 


5.3 


3.6 


2.8 


1.1 






K 


9.5 


6.9 


5.0 


2.4 


6 


8.0 


S 


2.3 


2.0 


1.2 


0.9 






P 


1.1 


1.0 


0.6 


0.5 






N 


1.4 


1.0 


0.7 


0.4 






K 


2.4 


1.9 


1.3 


0.8 




8.1 


S 


4.4 


3.7 


2.3 


1.7 






P 


2.1 


2.0 


1.1 


1.0 






N 


2.8 


2.0 


1.5 


0.7 






K 


4.7 


3.5 


2.5 


1.3 



a S=Salpetcr, F = 1.35; P=power law, T = 2; N=Naylor IMF, F = 0.2(M < 1M ) , F = 1.8(M > 1M ); K=Kroupa IMF,r = 0.3(0.080 > 
M > 0.5M Q ), T = 1.3(M > 0.5M Q ) 



